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Abstract 



A multivariate likelihood method to measure electroweak couplings with the Drell- 
Yan process at the LHC is presented. The process is described by the dilepton ra- 
pidity, invariant mass, and decay angle distributions. The decay angle ambigu- 
ity due to the unknown assignment of the scattered constituent quark and anti- 
quark to the two protons in a collision is resolved statistically using correlations 
between the observables. The method is applied to a sample of dimuon events 
from proton-proton collisions at ^/s = 7TeV collected by the CMS experiment at 
the LHC, corresponding to an integrated luminosity of 1.1 fb^^. From the dominant 
uu, dd — ?► 7*/Z — >• process, the effective weak mixing angle parameter is mea- 

sured to be sin^ 0eff = 0.2287 ± 0.0020 (stat.) ± 0.0025 (syst.) . This result is consistent 
with measurements from other processes, as expected within the standard model. 
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1 Introduction 

The main goal of the Large Hadron Collider (LHC) III is to explore physics at the TeV energy 
scale via proton-proton collisions. The Drell-Yan process fH occurs through the annihilation 
of a quark from one proton with an antiquark from the other proton, creating a virtual neu- 
tral gauge boson (7* or Z) that subsequently decays to a pair of oppositely charged leptons, as 
shown in Fig.[l] This annihilation process could reveal the existence of a new neutral gauge bo- 
son IHlIll or uncover deviations from the standard model of particle physics (SM) in elementary 
fermion couplings to the known neutral electro weak bosons. Besides quark-antiquark annihi- 
lation, new resonances decaying into lepton pairs may also be produced at the LHC via the 
gluon-fusion mechanism. 

Analysis of electron-positron annihilations into dilepton or bb pairs at LEP and SLC IHI led to 
high-precision measurements of the electroweak Z-boson couplings to fermions. The measure- 
ment of the weak mixing angle parameter sin^ 0w was performed to a precision of ~ 0.1%. In 
the SM, the weak mixing angle 0w describes the rotation of the original W*^ and B*^ vector boson 
states into the observed 7 or Z bosons as a result of spontaneous symmetry breaking [6|. Within 
the SM, sin^ 0w is the only free parameter that fixes the relative couplings of all fermions to the 
7 or Z. 

Measurements of the weak mixing angle with different initial and final fermion-antifermion 
states /i/i — 7- 7* /Z — > f2fi tests the universality of the fermion/ gauge-boson interactions and 
predictions of the SM. For example, LEP measurements [5J of the inclusive hadronic charge 
asymmetry provided a measurement of the couplings of light quarks compared to leptons and 
b quarks. The NuTeV collaboration measured sin^ 0w to precision of about 1% from neutrino 
and antineutrino deep inelastic scattering on nucleons [7J. The CDF and DO experiments |8]- 
[10 J at the Tevatron reached a similar precision with the Drell-Yan process in proton-antiproton 
collisions, as did the HI experiment IITll with electron-proton scattering at HERA. In this paper, 
we measure the sin^ Oy^ parameter to a precision ~ 1% in the proton-proton Drell-Yan process 
at the LHC with the CMS experiment. 

The proton-proton collisions at the LHC pose new challenges compared to previous collider 
experiments. The interference of the axial-vector and vector couplings leads to an asymmetry 
in the distribution of the polar angle of the lepton with respect to the direction of the constituent 
quark from the incoming proton in the quark-antiquark annihilation. This type of "forward- 
backward" asymmetry has been the primary measurement used to extract the couplings and 
sin'^ 0w at the LEP, SLC, and Tevatron experiments. However, because of the symmetric proton- 
proton collision at the LHC, the direction of the quark is not known and can be deduced only 
on a statistical basis ||4l [T2l - [T4l using the boost direction of the dilepton final state, because of the 
higher probability for a valence quark from one of the incoming protons to provide the boost. 
We have developed methods that allow a per-event analytical likelihood description to extract 
the maximal information about the process at the LHC. This technique exploits more informa- 
tion than the conventional forward-backward asymmetry approach, and the distribution of the 
polar angle as a function of both dilepton rapidity and mass is an essential component of this 
method. The technique has applications both for high-precision electroweak measurements 
and for rare-process searches at the LHC. 

In this paper, we present a multivariate analysis that uses the full information about the Drell- 
Yan process qq — > 7*/Z — ?> , parameterized as a function of the dilepton rapidity Y, 
invariant mass squared s, and decay angle 0* . This process offers a relatively simple environ- 
ment for the development of the matrix-element analysis techniques for resonance polarization 
studies at the LHC. Encouraged by feasibility studies of an analytical matrix-element approach 
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2 Phenomenology of the Drell-Yan Process at the LHC 



in Ref. flSl . we use a formalism based on an analytical description of the elementary interaction 
(Section |2]). The method is applied to a sample of proton-proton collisions at a center-of-mass 
energy a/s = 7TeV, corresponding to an integrated luminosity of 1.1 fb^^ recorded by the CMS 
experiment (Section |3]). We include a description of detector effects in the analytical likelihood 
model (Section |4]) and pay particular attention to systematic effects (Section |5|. The result is 
consistent with measurements in other processes, as expected within the SM (Section |6]). 





Figure 1: Left: diagram describing the SM process ^jij — > 7* /Z — > £ £+. Right: definition of the 
angle 9* in the production and decay of an intermediate state X, such as gg or (jij — > X — )■ i^£^. 



2 Phenomenology of the Drell-Yan Process at the LHC 

The philosophy of the multivariate likelihood analysis is to first produce a phenomenological 
model of the process and then introduce detector effects into the model. The parameters of 
the model may either be fixed to the best known values or left free in the fit, to be determined 
from data. These parameters may include the physical quantities of interest, such as sin^ 0W/ 
or a description of detector effects, such as a correction for the momentum scale in the track 
reconstruction. Therefore, we start with a discussion of a phenomenological model, and then 
proceed to detector-specific effects in the application of the analysis to CMS data. 

The tree-level coupling of a spin-one gauge boson to fermions is described by 

£^up^^e{pv- pA^s)vf, (1) 

where Vf and Uf are the Dirac spinors of the fermion (/) and antifermion (/), e'' is the polar- 
ization vector of the spin-one boson, and pv and Pa are the vector and axial-vector couplings. 
The couplings of the SM gauge bosons 7 and Z are given in Table [l] In the limit of negligible 
fermion masses, which is a good approximation for both quarks and leptons in the Drell-Yan 
process near the Z-boson mass, only two helicity states of the fermions are possible. They 
correspond to amplitudes oc [py — p^) and cx — {pv + Pa)- 

The parton-level cross section for the Drell-Yan process can be expressed with the help of the 
Wigner d^^^ ^, matrix, assuming the spin / = 1 intermediate states 7*, Z, and possible new 
unknown contributions (indicated by an ellipsis below), as 



%(s,cosr;0w) cx 1 X: (2/ + l)(C,A,-A,(^^* 



^x.x. ^x,M +^x.x. m^x.M (^w) X (T^^^^y^— ^ 



(2) 
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Table 1: Vector and axial- vector couplings of the SM gauge bosons to the charged fermion 
fields im. 



Pv Pa 



J —> e e+, ^ }i^,T T+ —1 

7^uu,cc,tt +2/3 

7^dd,ss,bb -1/3 

-3+12sin2 0w -1 



Z uu, cc, tt 
Z dd,ss,bb 



6sin(20w) 


2sin(2ew) 


+3-8sin2 0w 


+1 


6sin(29w) 


2sin(2ew) 


-3+4sin^0w 


-1 



6sin(20w) 2sin(2ew) 



where Xi'Xi, ^i, and A2 are the helicity states of the quark, antiquark, lepton, and antilepton, 
and the dilepton decay angle 6* is defined in the center-of-mass frame of the dilepton system 
as shown in Fig. [l] The effects of transverse motion of the incoming constituent quark and 
antiquark in their annihilation are minimized by using the CoUins-Soper frame lllTl . In this 
frame, the angle 6* is defined as the angle between the lepton momentum and a z'-axis that 
bisects the angle between the direction of one proton and the direction opposite that of the 
other proton in the dilepton rest frame. Should there be new states contributing to the cross 
section with the same or different spin, they will either contribute with interference if produced 
in qq annihilation, or without interference if produced in gluon fusion (gg) for example. 

The parton-level cross section in Eq. Q can be further simplified for the SM intermediate states 
7* and Z with spin / = 1 as 



d-qq {§, cos B*;ew) oc 



2(s^-m|)2 + m|r| 



3 [pr'Y (pV'^' 



2s 
2 



qq^Z 
Pv 



3(s-m|)pr'^V 
(s - m|)2 + m|r| 



X (1 +cos^r) 



qq^Z 
Pa 



pr" 



z^u 

Pa 



[1 + COS 



2 £1*^ 



pf^^pZ-^u (1 +cos2r) +2pf^^pr''cos0 



• (3) 



Equation (|3| implies that the asymmetry in the polar angle arises from terms linear in cos 9* . 
The value of the asymmetry depends on the dilepton invariant mass m = \/§, quark flavor 
q, and weak mixing angle 0w It could also be affected by deviations of couplings from SM 
expectations or by the presence of new contributions. 

The differential cross section of the proton-proton Drell-Yan process can be expressed as a 
product of the parton-level cross section in Eqs. (|2]) or (|3]l and the parton distribution functions 
(PDFs) fa{xi,s) [.18J describing the the probability for partons of type a to have a fraction Xj of 
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2 Phenomenology of the Drell-Yan Process at the LHC 



the proton momentum p = ^/s/2: 

dx,dx2dcose* " ? ((7,,(s,sgn(xi-X2)cos0;0w)/,(^i,s)/,(x2,s) 

(j=u,d,s,c,b ^ 

+ %(s,sgn(x2 - xi)cos0*;0w) fqi^lJ) fq{xi,s)^ . (4) 

The expression sgn(;ci — X2) refers to the sign of the difference {xi — x^), reflecting the fact that 
the quark direction is assumed to coincide with the boost of the qq system. This assumption 
introduces a dilution in the odd-power terms in cos d* . 

It is convenient to convert from the two variables {xi,X2) to {y,s), the dilepton rapidity and 
the square of the dilepton mass, as 

r = i.n(|±*Ul.nfilV (5, 



2 \E-pJ 2 \X2 
s = — p^ = X1X2S , (6) 

where E is the dilepton system energy, p and pz are its momentum and longitudinal momentum 
in the laboratory frame, and p = | Pz | at leading order in QCD. 

After transformation of the variables, the Drell-Yan process in proton-proton interactions can 
be expressed as follows: 

rfa-pp(Y,S,COSr;0w) T^even/. 9.*. X 

dYdsd cose* " ^ (s,cos ,M 



q=u,d,s,c,b 



+ D<j^(s, Y) X a°f^{s, cos 6*; 0w)J x F^^{s, Y) . (7) 

The parton factor is defined as 

y) = fq (e+ V^, i) fq (e-yiTs, s) + fq {e-WITs, s) (e+^v^, s) , (8) 
and the dilution factor is defined as 

f, [e+mV§rs, s) [e-m^/sTs, §) - fq [e-mVWs, s) fq [e+m^s, s 



The two components of the parton cross section contain either an even or odd power of cos 9* 
^7"(s, cos^ r ; 0w) = I (%(§, + cos 9*; 0w) + - cos 0*; 0w)) , (10) 



2 

a-°|^'i(s, cos 9*; 0w) = J + cos 9*; 9^) - %(s, - cos 9*; 9^)) . (11) 



The factors Fq^{s,Y) and Dg^(s,Y) both arise from the PDFs and can be extracted from pa- 
rameterizations such as in Refs. IIT9l - t22ll . We choose to extract PDFs numerically from the 
leading-order (LO) parameterization CTEQ6 ||T9]| to match the LO model of the process. We 
parameterize the PDFs analytically using polynomial and exponential functions in the relevant 
range of x and with coefficients that are also functions of s. The relevant range of x for this 
analysis is 1.1 x lO^'^ < x < 1.4 x 10^\ motivated by the detector acceptance and the dilepton 
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Figure 2: Analytical parameterization of the parton distribution functions xfa{x,Q^) at Q = 
100 GeV using the CTEQ6 [19J numerical computation for the various quarks, antiquarks, and 
the gluon. The gluon distribution is scaled by a factor of 0.1. 
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2 Phenomenology of the Drell-Yan Process at the LHC 




m (GeV) 

Figure 3: Distribution of the effective cross section factor 2mFqq{rrP- ,Y = 0) defined in Eq. (jsjl 
for five quark flavors (from top to bottom q = u,d, s, c, b) for proton-proton collision energies 
of 7TeV as a function of the dilepton mass m. An equivalent factor for gluon-fusion production 
is shown for comparison and is scaled by a factor of 0.1. 




Figure 4: The dilution factor D(j(j(s = m|,Y) for uu (top, red boxes) and dd (bottom, blue 
crosses) production as a function of the dilepton rapidity Y. The prediction from the PYTHIA 
simulation (boxes and crosses) of the 7*/Z — > ji^ji^ process and the analytical distribu- 

tions (solid curves) from Eq. ([9]l are shown. 




Figure 5: Distributions of Y (top), m (middle), and cos 6* (bottom), from Pythia simulation 
(points) of the —> 7*/Z — > process and analytical distributions from Eq. (|7jl. Distri- 

butions for five quark flavors are shown combined T,qq (black circles) and separately, in order 
of decreasing contribution: uu (red boxes), dd (blue crosses), ss (green diamonds), cc (magenta 
triangles), and bb (cyan stars). Distributions are normalized to unit area and are shown as 
fractions of events per bin. 
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3 Detection of the Drell-Yan Events with CIVIS 



mass selection criteria preserited in Section |4] This analytical parameterization is illustrated in 
Fig. |2] We choose an analytical parameterization because of the computational speed and ease 
of implementation. Its performance has been cross-checked with numerical computations. For 
systematic uncertainty studies, we check the performance against simulations with other PDF 
models 1111231, such as next-to-leading-order (NLO) CTEQ ||20l, MSTW (211, and nnpdf Il22l . 

The function F(^^(s, Y) is the effective cross section factor that scales the elementary parton-level 
cross section. This factor quickly decreases as the energy scale approaches values comparable 
to the full proton energies, as illustrated in Fig. Is] for Y = production. The factor D^^(s, Y) 
reflects the fact that the quark direction is generally unknown and is taken as the boost direction 
of the dilepton system, because of the higher probability for valence quarks to provide the 
boost. For q = u or d this factor ranges between and 1 as | Y| changes from to 4, as illustrated 
in Fig. |4] for s around the Z pole. From Eq. (|9| it follows that Dq^ = for q = s, c,b under the 
assumption fq{x, Q^) = fq{x, Q^), which is a good approximation in the current PDF model. 
A challenge at the LHC is that the dilution factor Dqq is small for the typical range of Y values 
in the detector acceptance region, as discussed below. Information about sin^ 0w or individual 
fermion couplings is contained in the shape of the three-dimensional distributions in Eq. ^ 
and enters through the elementary couplings of the electroweak bosons and fermions in the 
process ^ 7*/Z ^ 

Figure [s] illustrates the projections of the differential cross section from Eq. (|7]| in Y, m = \/s, 
and cos 9* for the five different quark flavors and combined. The relative fractions of the differ- 
ent quark flavors are 0.450, 0.375, 0.103, 0.048, and 0.025 for uu, dd, ss, cc, and bb, respectively. 
The results of the analytical model leading to Eq. ^ show good agreement with the predic- 
tions from conventional LO numerical Monte Carlo (MC) simulation using the PYTHIA genera- 
tor |24l with LO CTEQ6 1T9| PDFs. The above cross section is parameterized at leading order in 
both strong (QCD) and electroweak (EWK) interactions. Effects from NLO QCD contributions 
are studied with a detailed NLO POWHEG [25-27J simulation, which includes contributions 
from both initial-state gluon radiation and quark-gluon scattering. Effects from NLO EWK 
contributions are expected to be small compared to the precision of this analysis. EWK correc- 
tions are absorbed in a definition of the effective weak mixing parameter, sin^ 0eff [16]. We use 
sin^ ^eff in place of sin^ 0vv for the rest of this paper. 

We apply the above technique to the measurement of the weak mixing angle. We take the SM 
description of electroweak interactions and PDFs in the proton as given, and allow only the 
effective weak mixing angle 0eff to be unconstrained. More generally, the above formalism with 
the multivariate analysis of the Drell-Yan process allows us to study the elementary couplings 
of fermions to electroweak neutral fields, such as 7*/Z in the SM, as well as the structure 
functions of the proton. 

3 Detection of the Drell-Yan Events with CMS 

We apply the above method in an analysis of the qq^ 7*/Z — > process and measure 

sin^ 6f.fi. The choice of , as opposed to e^e+, is motivated by the more reliable descrip- 
tion of the detector and background effects, as well as the fact that this final state has not yet 
been studied for sin^ 0eff measurements in qq interactions. However, we do not expect any lim- 
itations in the method for future application to other final states. The expected multivariate 
distributions in Eq. ([7]| are modified by smearing due to detector resolution and photon final- 
state radiation (FSR), and by acceptance effects and non-uniform reconstruction as a function 
of the observables. All these effects are taken into account in the analytical parameterization. 
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as shown below. 

A detailed description of the CMS detector can be found in Ref. p8| . The central feature of 
CMS is a 3.8 T superconducting solenoid of 6 m internal diameter. Within the field volume 
are the silicon pixel and strip tracker, the crystal electromagnetic calorimeter (ECAL), and the 
brass /scintillator hadron calorimeter (HCAL). This analysis of the dimuon final state does not 
rely strongly on ECAL or HCAL measurements. Muons are measured in the window \f]\ < 2.5 
with the tracker and muon system. The pseudorapidity ?/ is defined as lncot(0/2) with the 
polar angle 6 measured in the laboratory frame. 

The silicon tracking detector (tracker) |29| consists of 1440 silicon pixel and 15 148 silicon strip 
detector modules. The pixel modules provide two-dimensional measurements of the hit posi- 
tion in the module planes, which translate into three-dimensional measurements in space, and 
are arranged in three layers in the barrel and two layers in the forward regions. The silicon 
strip detector is composed of 10 layers in the barrel region, four of which are double-sided, and 
12 layers in the endcap, where three out of six rings are with double-sided modules. Precise 
determination of the position of all silicon modules (alignment) is one of the critical aspects for 
achieving the designed resolutions of muon track parameters and is an important element of 
this analysis |30|. The muon system has detection planes composed of three distinct detector 
technologies installed outside the solenoid and embedded in the steel return yoke: drift tubes 
(in the barrel, \f]\ < 1.2), cathode strip chambers (in the endcaps, 0.9 < < 2.5), and resistive 
plate chambers (in both barrel and endcap regions, \t]\ < 1.6) |29|. 

This analysis uses data from proton-proton collisions at ^/s = 7 TeV collected during 2010 and 
2011, and corresponding to (1.07 ± 0.05 fb^^) of integrated luminosity. The signal and back- 
ground processes qq — > 7*/Z — > and t^t+ have been simulated with the NLO QCD 
generator POWHEG. Parton showering is simulated using Pythia. The NLO PDFs used are 
CTlO 1201. Background samples of W+jets and tt are generated using MadGraph [31 1, Pythia, 
and Tauola ||32|. Backgrounds from WW, WZ, ZZ, and QCD are generated using Pythia. 
Generated events are processed through the CMS detector simulation and reconstruction. The 
detector simulation is based on Geant4 Il33ll34l . 

Muon candidates are selected from a sample triggered online by events with at least two muons 
within the volume defined by |?/| < 2.4 and with transverse momentum (p^) requirements. 
These requirements depend on the period of data-taking; however, they always accept two 
muons with pj of at least 8 and 13 GeV, respectively. Offline, muon tracks are first reconstructed 
independently in the tracker and the muon system. Muon candidates are then reconstructed by 
two different algorithms [35] . The global muon algorithm matches tracks in the tracker to tracks 
in the muon system, and then refits the individual hits in the tracker and muon system to one 
overall track. The tracker muon algorithm extrapolates tracks in the tracker with > 0.5 GeV 
and p > 2.5 GeV to the muon system, and a track is taken to be a muon candidate if it matches 
at least one track segment in the muon system. Both algorithms take into account energy loss 
and multiple scattering in the steel yoke of the CMS magnet. Selection criteria demand at least 
10 hits in the tracker, including one in the pixel detector, at least one hit in the muon system, 
and a normalized < 10 for the global fit. 

Muons are required to have a small impact parameter, less than 2 mm measured with respect 
to the beam spot in the plane perpendicular to the beam direction. This requirement removes 
cosmic-ray muons and background events with displaced vertexes. We further require the an- 
gle between the two muon tracks to be larger than 2.5 mrad in the laboratory frame when the 
direction of one of the tracks is reversed. This removes any remaining cosmic-ray background 
and has negligible effect on the signal. To isolate single muons from muons overlapping with 
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4 Analysis Method 



jets, the sum of the transverse momentum of tracks in the tracker (excluding the muon in ques- 
tion) within a surrounding cone of AR = ^ (A?/)^ + {A(p)^ < 0.3 is required to be less than 15% 
of the measured transverse momentum of the muon, where At] and are the differences in 
pseudorapidity and in azimuthal angle in radians between the muon and the track. The ECAL 
and HCAL are not used for muon isolation, to reduce the effect from FSR and to maximize the 
amount of signal events. 

The kinematic requirements in the laboratory frame are \tj\ < 2.4 and pr > 18 and 8 GeV for 
the two muons. We introduce additional requirements in the Collins-Soper frame in order to 
simplify the acceptance parameterization: \tj*\ < 2.3 and pj > 18 GeV, where i]* and pj are 
defined with respect to the z'-axis, described previously. We also require the dimuon transverse 
momentum in the laboratory frame to be less than 25 GeV in order to suppress the contribution 
of events with hard jet radiation. Dilepton events are selected from events containing two 
oppositely charged, isolated, high-p^ muons with a dilepton invariant mass m in the range 
80 — 100 GeV. The dimuon rapidity Y is calculated from the lepton four-momenta as shown in 
Eq. (j5|. Restrictions on 6* and Y are motivated by detector acceptance effects, as discussed in 
SectionUl The number of selected events in the data is N = 297364. 



4 Analysis Method 

We use an unbinned extended maximum-likelihood fit that simultaneously describes the signal 
and background yields and the parameters of the (Y, s, cos0*) distributions. The likelihood 
function is written as 

£ = exp (-Msig - Mbkg) n ("sig X Psigixf, 0eff; + "bkg X 'Phkgixi) ^) j , (12) 

i 

where each event candidate i is characterized by a set of three observables Xj = { Y, s, cos 6*}i, 
Mgig is the number of signal events, which includes all intermediate states (7*, Z, and their 
interference), ftbkg is the small number of background events, Psigl^i;^) and V^kgi^i)^) are 
the probability density functions for signal and background processes, and ^ represent the 
parameters of these functions. The signal probability density function is defined as 

/+00 
dx n{x) Videa\{y, S - X, COS 6*; 0eff) • (13) 



The ideal distribution VideaiiY, s, cos 6*; 0eff) in Eq. ( 13 1 is the Drell-Yan cross section defined in 
Eq. 0. We correct 'Pideal for detector effects, such as acceptance, parameterized with Q{Y,s, cos 9*), 
and resolution and photon emission (FSR), parameterized with TZ{x) where x is the change in 
the dilepton center-of-mass energy squared. 

The acceptance function ^(Y, s, cos0*) describes the nonuniform reconstruction efficiency as 
a function of the three observables, which includes effects from online trigger requirements, 
detector acceptance, reconstruction algorithms, and selection requirements. The most impor- 
tant effect is the loss of particles near the beam directions and the second-most-important effect 
is the minimum transverse momentum requirement on the leptons. Otherwise, the efficiency 
across the acceptance range, defined by the selection requirements \?]*\ < Y^ax = 2.3 and 
Pt > Pmin = 18 GeV, is close to uniform. The above selection requirements define a sharp 
boundary in (Y, s, cos0*) space, which can be expressed as limits on cos 6* for given Y and s 
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values as follows: 



I cos 0*1 < tanh(y„ 
I cos( 



1^1). 



This boundary is illustrated in Fig. |6]in the [Y, cos 9* ) plane for a fixed value s 
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Figure 6: Accepted cos 6* range as a function of Y for s 
used in this analysis. The outer boundary corresponds to Eq. ( 14 1 and the horizontal lines near 
COS0* = ±0.92 correspond to Eq. (jlSl). 



The effect of smearing the muon track parameters, such as the muon momentum and direc- 
tion, due to detector resolution and FSR, is most evident in the invariant mass distribution. 
This effect is parameterized with the function TZ{x) in Eq. ( |T3) |. Both acceptance and reso- 
lution effects are illustrated in Fig. [zj where the analytical parameterization of Y, m = \/s, 
and cos 6* is in good agreement with LO simulation in both QCD and EWK, as generated 
by PYTHIA. Although a wider m range is investigated, the analysis is performed in the range 
80 < m < 100 GeV to reduce uncertainties from FSR. In this illustration, FSR is included and the 
major detector effects are introduced in the following way: for the three track parameters (pr, (p, 
9), we apply Gaussian random smearing with standard deviation of Apj = O.OlBpr + 0.0001 
(with in GeV), = A9 = 0.001 rad, and neglect resolution effects on the track origin. This 
simplified simulation of detector effects is found useful to isolate production model uncertain- 
ties from the detector effect parameterization. 

Further studies are performed with full GEANT4-based modeling of the CMS detector using the 
POWHEG simulation of the dimuon events and with PYTHIA simulation of the parton show- 
ering and FSR. In the parameterization of the acceptance function Q{Y,s,cos9*), we model 
the small deviations from a uniform efficiency with empirical polynomial functions that in- 
clude correlations of the two observables within the boundaries of the (Y, cos 9* ) plane defined 
above. This efficiency parameterization is derived from the simulation with a fit where the 



12 



4 Analysis Method 
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Figure 7: Distributions of Y (top), m (middle), and cos 9* (bottom), from Pythia simulation 
(points) of the qq 7*/Z H process and its analytical parameterization (smooth 
curve). Combined distributions from Fig. |5] appear in red at the top of each plot (red circles 
for an "ideal" simulation), while distributions after acceptance and resolution effects, includ- 
ing photon FSR, appear in blue below (blue squares for a simplified "detector" simulation). 
Distributions are normalized to unit area and are shown as fractions of events per bin. 
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parameters of the polynomial functions are left unconstrained. The main effect is a loss of ef- 
ficiency in the vicinity of the acceptance boundaries. A similar approach is later employed as 
part of the systematic uncertainty studies where the parameters of the efficiency model are left 
free in the fit to data. 

In the parameterization of the resolution function TZ{x), FSR is modeled with Pythia and res- 
olution effects are taken from the full CMS detector simulation, including the effects of tracker 
alignment on the tracking resolution. The function TZ{x) is approximated with a sum of four 
Gaussian functions, to allow for the analytical convolution in Eq. p3] l and be flexible enough to 
describe both detector resolution and FSR effects. Parameters of the Tl{x) function are left free 
in the fit to the simulated MC sample. The overall shift of the Z mass in the resolution function 
TZ{x) is left free in the fit to data, effectively allowing the energy scale to be determined from 
the data. 

The background contribution is estimated by MC simulation; the QCD component has been 
cross-checked with data. The total expected background is about 0.05% of the signal yield. 
The background consists of the crossfeed from the qq Z/7* t+t^ process, QCD, tt, 
and diboson production in nearly equal contributions. The probability density function for the 



background Vhkgi^i', ^) is parameterized in a similar manner to that shown in Eq. (13i with 



an acceptance range defined by Eqs. (|T4)| and ( 15 1 and the distributions within the acceptance 



boundaries parameterized with an empirical polynomial function. The number of background 
events n^kg is fixed to the expected value of 157 events. 

In Fig. [s] we show the cos 9* distributions in the data separately for the |Y| < 1 and |y| > 1 
regions, and compare them to the POWHEG-based simulation of the qq^j*/Z^}i^}i^ pro- 
cess in the CMS detector. Together with Fig. |4| these distributions illustrate the challenge of 
analyzing Drell-Yan events at the LHC. While the acceptance effects on the cos 6* distribution 
are moderate for smaller values of | Y|, the dilution is strong, as shown in Fig. [Ij In contrast, 
the larger values of | Y | have a smaller dilution effect, but the cos 9* range is strongly truncated 
because of the limited acceptance, as shown in Fig. |6] Therefore, the optimal analysis of the 
angular distributions requires proper accounting for such correlations among the three observ- 
ables. At the same time. Fig. |8] shows good agreement between the data and MC simulation. 
Residual differences in the distributions can be explained by the somewhat different value of 
sin^ ^eff = 0.2311 used in the simulation compared to the best value describing the data, and by 
several systematic effects accounted for below, such as the tracker misalignment, the momen- 
tum scale in the track reconstruction, and FSR modeling. 

A "blind" analysis of the data is performed, in which the fit result is not examined until a review 
of the entire analysis is complete, including the evaluation of all associated systematic uncer- 
tainties. However, while the analysis is performed "blind," the quality of the fits to the MC 
simulation and data is examined. We test the performance of the fitting procedure using sam- 
ples generated using Monte Carlo simulations, with each separate sample containing the same 
number of events observed in the data. Signal events are generated with the PowHEG-based 
CMS detector simulation with an input value of sin^ 0eff = 0.2311. The number of background 
events is Poisson distributed according to expectation. After the corrections discussed below 
are applied, the pull distribution is in agreement with a unit-width Gaussian distribution cen- 
tered at zero. A comparison of the MC sample projections and the probability density functions 
are shown in Fig. |9] 

We examine the quality of the fit to the data by comparing the data distributions to the likeli- 
hood model expectations, and comparing the fit likelihood value C and the observed statistical 
uncertainty to those expected with the generated samples. Projections of the data and the prob- 
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ability density functions are shown in Fig. 10 They exhibit similar agreement as with the simu- 
lation shown in Fig. |9] Correction for the energy scale is already included in the fit model. The 
observed statistical error on sin^ 0eff of ±0.0020 is in good agreement with what is expected 
from the MC samples discussed above. We find only small differences when comparing the 
likelihood value L for generated experiments from the likelihood model, the PowHEG-based 
CMS detector simulation, and the data. The level of agreement is consistent with typical dif- 
ferences due to imperfect efficiency function modeling and NLO effects discussed below. The 
variations do not affect the result of the analysis within the systematic uncertainties assigned. 



5 Systematic Uncertainties 

The list of systematic uncertainties on the measurement of sin^ 0eff arid associated corrections 
to the fit values, as discussed below, is shown in Table |2) These uncertainties arise from both 
theoretical assumptions and detector modeling. 

Table 2: Corrections to the fit values and systematic uncertainties in the measurement of 
sin^ 0eff ■ 



source 


correction 


uncertainty 


PDF 




±0.0013 


FSR 




±0.0011 


LO model (EWK) 




±0.0002 


LO model (QCD) 


+0.0012 


±0.0012 


resolution and alignment 


+0.0007 


±0.0013 


efficiency and acceptance 




±0.0003 


background 




±0.0001 


total 


+0.0019 


±0.0025 



We follow the PDF4LHC working group's recommendation | |23| in estimating uncertainties 
from the PDFs. We reweight a large MC sample generated with CTlO [20J PDFs to obtain 
samples equivalent to MSTW 2008 [21] and NNPDF 2.1 [22J PDFs. We vary the internal degrees 
of freedom of the PDFs for all three sets of models. We also use this technique to vary dls, but 
find its uncertainties to have negligible effects compared to any of the PDF variations. We find 
a change in the value of sin^ 0eff of +4-8 x 10^^ and ±3.4 x 10^^ in using the PDFs from MSTW 
and NNPDF, respectively. The variations in the value of sin^ Sgff within each MC simulation due 
to the PDF uncertainties are x 10""^, x 10""*, and ±7.3 x lO^'^ for CTlO, MSTW, and 
NNPDF, respectively. The envelope of the above variations corresponds to the total systematic 
uncertainty of ±0.0013. 

The FSR is modeled with Pythia in the parameterization of the resolution function IZix). 
As a cross-check, we use four alternative FSR models for generation, with simplified detector 
simulation discussed above: PYTHIA, Photos [36], and two different modes in HORACE [l37|. 
All three generator programs perform 0(a) calculations of FSR and provide similar results, 
leading to differences in the fitted values of sin^ Sgff of about 0.001. In addition, the HORACE 
generator allows the exact 0(a) calculation and multiple-photon radiation from all charged 
states, which is the state-of-the-art EWK calculation. We found that this has a larger effect on 
the analysis when a wide range of values for the dimuon invariant mass m is used. However, 
with the mass range 80 < m < 100 GeV, the differences in the relevant part of the radiative tail 
become small. We perform cross-fits of the four generated samples and the four correspond- 
ing resolution functions 7?.(x), finding differences in the fitted sin'^ 0eff values of at most 0.0011. 
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Figure 8: Distribution of cos 6* for data (points) and PowHEG-based detector simulation with 
sin^ 0w = 0.2311 (histogram) of the qq —)■ j* /Z —> process for |y| < 1 (top) and |y| > 1 

(bottom). 
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MC sample corresponds to an integrated luminosity of 16 fb^^. The lines show the projections 
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The Pythia sample typically results in larger differences from the other generators, and the 
assigned systematic uncertainty of ±0.0011 covers these deviations. The assigned error conser- 
vatively covers the FSR uncertainty in the measurement when the fraction of FSR radiation is 
reweighted in simulation. This reweighting technique is based on the comparison of the FSR 
fractions between data and the Pythia simulation. 

Effects from NLO EWK contributions are expected to be small compared to the statistical pre- 
cision of our measurement. Calculations with the Zfitter 1138 1 program indicate that the vari- 
ation of the effective value of sin'^ ^gff for light quarks is within 0.0002 of the lepton values. It 
is only the heavier b quark that requires corrections of the order 0.001. However, given that 
only about 2.5% of the dileptons are produced in bb interactions, and no angular asymmetry 
can be measured because of the dilution factor with this initial state, these corrections have a 
negligible effect on our measurements. 

Parameterization of the likelihood function models both the initial-state interactions and the 
PDFs at LO. The requirement that the dimuon transverse momentum pj be less than 25 GeV 
suppresses the contribution of events with hard jet radiation and reduces the effects from NLO 
processes. This requirement also ensures that the transformation between the laboratory frame 
and the Collins-Soper frame is small, and that our analytical description of acceptance effects 
is correct, without any loss of acceptance coverage. With generated samples, we observe a 
bias of -0.0012 ± 0.0006 in the fit value of sin^ e^n, which is attributed to NLO effects. In this 
test, perfect CMS detector conditions are simulated, which removes most of the detector effects 
discussed below. To be conservative, we apply a correction of +0.0012 and assign a systematic 
uncertainty of ±0.0012 to cover all effects associated with the LO model. For example, we have 
investigated the dependence of the expected and observed shift in the fit result as a function of 
the Pt requirement over a wide range of pr- The results are stable within the uncertainty quoted. 
The distribution of pr exhibits only a small difference between the data and the MC simulation. 
We reweight the MC simulation to match the pj distribution to the data and observe that the 
results of the fit to the reweighted MC events are consistent with the case without reweighting 
to within 0.0004. Treating the correction as an additive or multiplicative factor does not affect 
the final result, as long as the observed value of sin^ Sgff is close to what is expected. 

The detector resolution in the muon track reconstruction is affected by contributions from the 
silicon tracker alignment. We perform a realistic simulation of the alignment procedure to 
model the statistical precision of the track reconstruction. We observe a bias of —0.0013 in the 
fit result of sin^ 0eff when the realistic simulation of alignment is used in place of perfect con- 
ditions. We find that while the statistical precision of the track reconstruction is well modeled 
by the realistic simulation, the biases from -invariant detector deformations [30J may not 
necessarily be well reproduced in the MC simulation. We have investigated nine basic distor- 
tions in the tracker reconstruction geometry, which follow from the cylindrical symmetry of 
the system ||30ll39l. In each case, the procedure of the tracker alignment is repeated after the 
distortion is introduced. The effects of the remaining distortions on sin^ 0eff ^re all smaller than 
0.001. The typical initial distortions are taken to be 200 f/m, which is the approximate value of 
the constraints from the detector survey, the placement tolerance, and the observed agreement 
in the alignment procedure. 

In the end, constraints on the above distortions in the tracker reconstruction geometry come 
from the data. We have observed that distortions affecting the sin^ 0gff fit values also introduce 
a bias in the mass of the dimuon pair m as a function of cos 6* . We observe a linear trend in 
the bias of the average value of m as a function of cos0*, with a slope of —0.072 GeV, in the 
realistic simulation that is twice as large as that observed in data, —0.039 GeV, when both are 
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compared to simulation with the ideal geometry model. We also observe a bias in the value 
of sin'^ 0eff that is twice as large when an additional systematic distortion is introduced in the 
realistic simulation, resulting in the slope of the average m value versus cos d* dependence also 
becoming twice as large. From these studies, we assign a correction of +0.0007 to the fit value 
of sin^ 0eff due to alignment effects and a systematic uncertainty of ±0.0013 to cover the range 
of possible deviations observed. In order to minimize the uncertainties from the energy scale 
bias in the track reconstruction, the shift of the Z mass in the resolution function TZ{x) is left 
free in the fit, effectively allowing the energy scale to be determined from the fit to the data. 
Consistency between the fit value from the data and the expectation from the MC simulation is 
found to be within 0.1 GeV. 

We find very weak sensitivity to the efficiency parameterization Q ( Y, s, cos 9* ) across the ac- 
ceptance range because the efficiency is symmetric in cos 0* . This leads to negligible effects 
on the odd terms in the angular distribution that are sensitive to sin^ Sgff- The sign of Y is de- 
fined by the dimuon system direction along the counterclockwise beam and has no preferred 
direction. The sign of cos 9* is defined by the charge of the "forward" lepton. The cylindri- 
cal symmetry of CMS, combined with the random nature of the "forward" direction, leads to 
a symmetry in the efficiency function. This has been verified with a detailed GEANT4-based 
simulation of the CMS detector, including calibration and alignment effects. Even in the ex- 
treme case of ^(y,s,cos0*) being flat across the acceptance range, negligible changes in the fit 
results are observed with simulated samples. We also allow parameters of the model to be free 
in the fit to data. We assign a systematic uncertainty of 0.0003 due to efficiency and acceptance 
parameterization, which is the level of consistency of results from these studies. 

The number of background events ntkg is fixed to the expected value and is varied according 
to its associated uncertainties. We assign a 50% uncertainty to the QCD rate, based on studies 
with wrong-sign lepton pairs. The relative size of the sum of the EWK background processes 
is expected to be reproduced by simulation to a precision of better than 20%. However, in the 
mass range 80 < m < 100 GeV, the fraction of background is only 0.05%, and the fit results 
are insensitive to the exact treatment of the background. The measured sin^ ^gff value remains 
stable within 0.0001, even when the background is removed from the model. 

6 Results and Discussion 

We have presented a likelihood method to analyze the Drell-Yan process at the LHC. The pro- 
cess is described by the correlated dilepton rapidity, invariant mass, and decay angle distribu- 
tions. The quark direction in the elementary parton collisions, which is not directly accessible in 
the proton-proton collisions at the LHC, is modeled statistically using correlations between the 
observables. The result of the analysis, which includes systematic uncertainties and corrections 
from Table |2} is 

sin^ = 0.2287 ± 0.0020 (stat.) ± 0.0025 (syst.) . 

This measurement of the effective weak mixing angle in the predominantly uu, dd — ?► 7* /Z — ?► 
processes in proton-proton collisions is consistent with measurements in other processes 
MSllTHTTl, as expected within the standard model. 

The dominant systematic uncertainties in the measurement include modeling of the PDFs, FSR, 
effects beyond the leading order in QCD, as well as detector uncertainties primarily due to 
tracker alignment. With increased statistics of the Drell-Yan process at the LHC, a further 
reduction of the systematic uncertainties will become critical. Understanding the tracker align- 
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ment will certainly improve as the collaboration gains further experience. Therefore, we expect 
the limiting imcertainties to come from the Drell-Yan process modeling. 

Uncertainties from PDFs will decrease as better constraints on the proton model become avail- 
able from the LHC and elsewhere. In fact, the Drell-Yan process is itself a useful input to the 
PDF model constraints, and the methods discussed in this paper can be used to constrain the 
parameters in the PDF model. However, one must be careful not to mix information used for 
PDF constraints from the Drell-Yan process with measurements using the same events, unless 
the correlations are properly taken into account. Uncertainties from the FSR model may be im- 
proved as higher-order electroweak calculations are integrated with the higher-order QCD cal- 
culations of the matrix element in the Drell-Yan process, such as the incorporation of POWHEG 
and Horace. 

The LO approximations in the model may be further improved as NLO matrix elements are 
employed in the likelihood approach and more variables are integrated into the analysis. We 
view the current LO formalism as a conceptual step in developing multivariate matrix-element 
approaches to resonance polarization analyses, which can be applied to precision measure- 
ments, as well as potential new resonances that may be discovered at the LHC. The evolution 
of this method may also allow several parameters of the electroweak couplings to be deter- 
mined simultaneously, such as a measurement of the vector and axial-vector couplings of the 
light quarks separately from the lepton couplings. 
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